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Abstract. - One or more small holes provide non-destructive windows to observe corresponding 
closed systems, for example by measuring long time escape rates of particles as a function of hole 
sizes and positions. To leading order the escape rate of chaotic systems is proportional to the hole 
size and independent of position. Here we give exact formulas for the subsequent terms, as sums 
of correlation functions; these depend on hole size and position, hence yield information on the 
closed system dynamics. Conversely, the theory can be readily applied to experimental design, for 
example to control escape rates. 



A fundamental issue in physical measurement (whether 
classical or quantum) is ensuring that the system to be 
studied is little affected by the observation. This is clearly 
much easier to satisfy if the measuring device lies outside 
, the system. A truly isolated system cannot affect its sur- 
roundings, and is impossible to observe. However it can 
be to a good approximation unaffected by a very small 
hole, through which particles or radiation can carry infor- 
mation about the internal dynamics. Notice that theory 
and experiment often work in opposite directions: theo- 
rists use closed systems to understand open ones, while 
experimentalists do the reverse. 

The initial motivation for this work came from atom- 
optics billiards [1,2] in which escape properties were used 
to distinguish regular and chaotic behavior, but for which 
a detailed theory was absent, a deficit which the current 
work aims to address. Our theory is however far more gen- 
eral than these experiments suggest; billiards in which a 
particle moves in straight lines making specular reflections 
with a boundary are relevant to any system with particles 
in a relatively homogeneous cavity. Experimental real- 
izations have included microwaves [3-6], visible light [7], 
phonons in quartz blocks [8] , and electrons in semiconduc- 
tors [9,10]. Furthermore, transport through small holes (in 
phase space) has applications as diverse as the transition 
state theory of chemical reactions [11], the migration of 
asteroids [12], and passive advection in fluids [13]. It 
can also help to characterise chaos in relativity [14] and in 



Hamiltonian systems in general [15]. Finally, the "hole" 
can be a desired region in phase space, for example a set 
from which we can subsequently control the dynamics with 
small perturbations; the escape rate then gives informa- 
tion on how long we may need to wait before attempting 
control [16]. We emphasize that while our numerical ex- 
ample is a billiard, the method is readily applicable to 
general volume preserving (for example Hamiltonian) sys- 
tems. 

Just as black body radiation emanating from a cavity 
through a small hole gives only a single piece of infor- 
mation (the temperature), we would expect the transport 
properties of orbits in the phase space of open systems 
to give only simple geometrical information to leading or- 
der. The transit time (for orbits entering through a hole 
then exiting) is exactly given by the ratio of hole to sys- 
tem sizes [17]. The escape rate (for orbits initially in the 
system) is hard to characterize in general [18], but is also 
often assumed to follow a similar equation [17]. In a pre- 
vious Letter [19], we showed a result of this type for the 
regular circle billiard, with corrections computed analyt- 
ically. Stadium and related billiards with interior 
holes have also been studied recently [20]. Here we 
consider general strongly chaotic systems in the limit 
of small holes, and find an analytic expansion of the es- 
cape rate, in which the ratio of hole and system sizes gives 
the leading term, and corrections are given by correlation 
functions. Measurement of the escape rate of a hole as 
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a function of size and position, or of two holes compared 
to the individual escape rates, then provides a window 
through which to study the original (closed) dynamics. 

The use of correlation functions to compute transport 
coefficients is familiar from Green-Kubo formulas [21]. 
Also, the escape rate from rather specific systems of large 
spatial extent can be used to compute general transport 
coefficients [22] . Thus the appearance of correlation func- 
tions in more general escape rate calculations should not 
come as a surprise. We also note that the periodic orbit 
formalism for computation of escape rates of open chaotic 
systems [23,24] is unsuitable for the limit of small holes, as 
a very large number of periodic orbits would be required 
for the hole(s) to be sufficiently covered. 

Our formalism is in the setting of volume preserving 
maps $ : M — » M, such as Hamiltonian evolution con- 
sidered at equally spaced times (a "stroboscopic map" ) or 
when a certain condition is met (a "Poincare map"). We 
use (/) to denote an integral of a phase variable / : M — > M 
over the normalised volume clement of the phase space M. 
Even when the dynamics does not preserve the 
usual volume (for example leading to a strange at- 
tractor), there is often a well defined fractal mea- 
sure which is preserved, albeit making calculations 
of correlation functions more complicated. 

A particular case of a Poincare map is that of a billiard 
system of a particle moving uniformly between specular 
collisions with the boundary 3D of a two dimensional do- 
main D, where $ denotes the evolution from one collision 
to the next. In this case the phase space volume element 
of the flow projects to the area element dldp\\ of M pre- 
served under $, where I measures arc length along the 
boundary and p|| is the component of momentum paral- 
lel to the boundary after the collision. Thus for a phase 
function / on the billiard boundary we have 



J f(l, Pll )dldp l{ 
J dldp\\ 



(1) 



where the denominator is simply 2|p||9£)|. 

We now return to the general formalism and define some 
phase variables on M . Let the function T : M — > K denote 
the time from one collision to the next; this is necessary to 
relate collision number with real time t. We assume this 
function is bounded from above. The holes are represented 
by a characteristic function \ '■ M — > M equal to zero on 
a hole and one elsewhere; in the billiard example this is 
normally a function only of I and not p\\ . Note that our 
approach here is to always use the dynamics of the closed 
system, but "kill" escaping trajectories using a multipli- 
cation by x- (x) = 1 — h is simply the fraction of M not 
covered by holes and h is the size of the hole (relative to 
|M|). Another useful set of averages is r k = (T k \) / (x), 
the average of the fcth power of the collision time over 
the non-hole part of the system. We will need a weighted 
characteristic function Xs '■ M — > M defined by Xs = e sT x, 
so that xo = X- It is also useful to define phase functions 



with averages subtracted; we find that the best way to do 
this is to set Xs = Xs/ (Xs) - 1 so that x = Xo = x/(x> - 1- 
We also define f = T - n and r k = (f k x )/(x)- 

Assume that the phase space is filled with initial con- 
ditions (or mutually non-interacting particles) of a uni- 
form density with respect to the preserved volume ele- 
ment, then consider the survival probability P{t) for time 
t, so -P(O) = 1. The exponential escape rate 7, defined 
for strongly chaotic systems, is 



7 = — lim 

t — >oc 



lnP(t) 
t 



(2) 



which is well behaved in strongly chaotic systems and is 
equal to the leading pole (that is, of smallest s > 0) of the 
function 



P(t)dt 



f 
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e st (xV ...X Nt )dt 



(3) 



where N t is the number of collisions before time t and a 
superscript on a phase variable from now on will usually 
denote discrete time, ie x 3 = X $ 1 1 1 with the map 
$ composed j times. We can rewrite this as 



f 

Jo 



N=0 



where Xs is defined above, and we have ignored a factor 
bounded in magnitude which does not affect the leading 
pole. Now assuming that the long time escape is not dom- 
inated by orbits with very small T (hence T N is also ef- 
fectively bounded), the escape rate is given by the leading 
pole of 

00 00 
G( S )=^G N ( S )=^(x°x^...xf) (5) 
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Substituting Xs = (Xs)(l + Xs) to isolate the leading part, 
and then expanding and grouping into terms with different 
numbers of Xs, we find 



G( s ) = ^(x s ) w+1 r N ( s ) 



(6) 
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where 



N+l 



= E E (si 1 *? ■■■«"> (7 

n=0 0<j 1 <j 2 ...<j n <N 

= 1+ E (xix k s }+ E (xixsxi> + -- 

0<j<k<N 0<j<k<l<N 

Remarks: 

1. The single sum (n = 1) is absent since the average of 
Xs is zero by definition. 

2. This expansion is divergent since the number of terms 
in each sum gets larger with each term, and the higher 
order correlations are small but finite. 



Peeping at chaos 



The sum for G(s) diverges when the ratio of subsequent 
terms exceeds unity in the limit, so at the first zero of 



9(s) 



lim g N (s) 

N—>oo 



(8) 



where 



g N (s) = In 



G 



N+l 



(s) 



ln( Xs )+lnr^ +1 (s)-lnr JV (s) (9) 



G N (s) 

Expanding the logarithms and taking the limit, we find 



n=2 



g{s) = ln( Xs ) + Q™( s ) 
where the cumulants Q n {s) are 

Q2{S) 

Q 3 (s) 



(10) 



^ ^ (XsXs) 
^ ^ (XsXsXs) 



(11) 



o<j<fc 



The higher cumulants are more complicated, for example 
Qi(s) includes 4-time correlations and products of 2-time 
correlations. We expect this cumulant expansion (by anal- 
ogy with [25]) to be well defined if the system has multi- 
ple correlation functions decaying faster than any power 
(which, a posteriori, is what we mean by "chaotic" ) . De- 
cay with a power greater than unity will permit Q2(s) to 
exist, and hence allow the second order formulas below to 
be used. 

We are interested in the limit of small holes, ie small 
escape rate. Thus we expand these quantities in powers 
of s: 



Xs 



x + (i + x) 



sT- 



(*•-*) 



(12) 



Note that the first term is independent of s and the re- 
mainder is projected onto the non-hole part by 1 + \- 



ln(x s ) = ln(x) + sti + —t 2 + . 



(13) 



Note that higher terms in these expansions involve prod- 
ucts of the various fk , and for Xs , powers of T as well. The 
cumulants can similarly be expanded, but it turns out we 
will need only the leading ,s = term, that is, replacing 
Xs by x- 

We need to establish the order of magnitude of the var- 
ious quantities in the expansion. First we note that ln(x) 
is of order h (the hole size). The powers of T and their 
correlations are of order unity. Each factor of x in a cor- 
relation will result in a factor of h unless there is a high 
probability for orbits entering the hole to return there. 
In the case of a billiard in which a short periodic orbit 
lies on the hole, an orbit leaving the hole needs to be in 
a precisely specified direction in order to return once in 
a small number of collisions, but after this it will return 



with high conditional probability. Thus correlations with 
two or more x terms should be of order h 2 . In non-generic 
billiards, (for example if the boundary contains an arc of 
a circle centred on the hole that reflects orbits in many 
directions back to it) or for some non-billiard systems, no 
fixing of direction is necessary and all correlations could 
be of order h. 

More generally, the first order at which periodic orbits 
contribute depends on the dimension of the system, the 
size of the hole in each dimension and the structure of 
the periodic orbits (eg whether isolated). Typical low di- 
mensional chaotic systems have periodic orbits that are 
dense, so all holes will cover sufficiently long periodic or- 
bits. However the probability of following a long orbit for 
a whole period is very small, related to the exponential of 
the Lyapunov exponent times the period. Hence the term 
"short periodic orbit" means that the set of orbits close to 
them have a probability of return large enough to make a 
measurable impact on the escape rate. Previous work [16] 
has treated holes on short periodic orbits, but limited to 
first order in h. 

We note immediately that if we can assume the cor- 
relations of x are small, the escape rate 7 given by the 
solution of g(s) = reduces simply to 7 = — ln( X )/ri, 
which is the hole size divided by the average collision time 
to leading order, as expected. It also means that s is of 
order h. Now, given our assumptions, we know the order 
of each of the quantities and we can proceed to develop 
an expansion in the single parameter h. 

We expand \s into orders in h: \s = Xs^ + X^ + • • • 
with 



<P = X + sf(l+x) 



= y(T 2 -i 2 )(l + x) 



(14) 



Note that while sTx in X^ 1 ^ appears to be of higher order, 
it needs to be included so that (xi^) = at each order; 
this is essential for convergence of the correlation sums. 
g{s) now splits into terms based on their order in h, taking 
care that in general, x behaves differently to s: q(s) = 
g( 1 \s)+gW(a) + ... with 



g (1 \s) = ln(x)+sn 



(15) 



g (2) (s) = ^f 2 + ^(xW°xi 1)J ) + g3(o) + g 4 (o) + ... 

0<j 

Higher order contributions arc more complicated, and 
most significantly, all involve an infinite series of corre- 
lations. Far from a short periodic orbit, we can hope 
that /c-order correlations are of order h k , leading to a 
more tractable expansion, with higher order cumulants 
relegated to higher orders in s. 

Finally the solution of g(s) — gives the escape rate 7. 
Writing (in orders of h) 7 = 7W +7^ + . . ., expand 3(7) 
in a Taylor series about 7^. Thus we have 



v(!) 



ln(x) 



„( 2 ) 



l(2) 



( 7 (D) 



(16) 
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This confirms the leading order behavior that the escape 
rate is the hole size divided by the average collision time, 
and observe that the result involves correlations of a quan- 
tity 



y (1) 



x- 



ln( X > 



T(l + X ) 



which to leading order is equivalent to 



h 



(17) 



(18) 



Thus we arrive at our first main result, giving an exact 
expansion for the escape rate in powers of the hole size, 



1 

n 



ln(l -h) + ^ 



ln(l - h) 



n 



n—3 



+ 5>V) 

3=0 

(19) 



where the higher cumulant terms are ignored for the pur- 
poses of the simulations below and in general if the hole(s) 
avoid short periodic orbits. If the hole(s) lie on short 
periodic orbits, relevant terms in the higher cumu- 
lants would need to be summed explicitly. This is 
beyond the scope of this Letter, but should repro- 
duce and generalise relevant equations in Ref. [16]. 

Let us now compare the escape rates of a billiard with 
one hole of size a and corresponding X A, one (disjoint) 
hole of size b and corresponding \b, and both holes, total 
size a + b and corresponding xab = Xa + Xb — 1- We 
compute (up to second order in the hole size) 



(2) 

Tab 



Ji) , ri) , ab HTxa) _ a{T XB ) 

I A ' IB ' lr^\ /rr,\ O 1^1 O T 



(2) , (2) 

Ta +Tb 



(T) (T)2 
abi 2 



(T)3 



(20) 



0<j 



+ iQnAB(0) - Q nA (0) - Q„fl(0)] \ + ... 
n=3 J 

where u for a single hole {A or B) is defined above. Note 
that no subscript need be given for T or f 2 since to leading 
order they are equivalent to the ones of the closed system, 
namely f = T - (T) , f 2 = (T 2 ) . Putting this all together 
we find a remarkably simple relation, our second main 
result 

1 f °° 

1AB = 1A + IB - jp- < E ("A«B> 
[j=-oo 

oo ^ 

+ E iQnAB(O) - Q nA (0) - Q nB (0)} }+■■■ (21) 




n=3 



r = o 

Fig. 1: The billiard with R = l/y/3 showing holes of size R/20 
at positions used in Fig. 2, and the I* coordinate used in the 
later figures. 



The higher correlations all contain mixtures of xa and xb, 
so we expect they contribute at second order only if both 
holes lie on the same short periodic orbit. In Ref. [16] 
the difference jab — 7a — 7s is termed the "inter- 
action" between the holes, and described approx- 
imately in terms of "shadows" cast by one hole 
on another, that is, overlap between one hole and 
the image of another under the dynamics. Here 
we give a precise formulation in terms of correla- 
tion functions, which could also be extended to the 
three and multi-hole interactions described there. 
Clearly the strongest interactions exist when there 
are large shadow effects at short times. Whether 
this is the case or not, the full (long time) corre- 
lation functions permit a precise calculation. 

Finally, we test the above formulas with numerical sim- 
ulations. We consider a "diamond" billiard, bounded by 
four overlapping disks of radius R, centered at the cor- 
ners of the unit square (Fig. 1). When 1/2 < R < l/y/2 
this has strong chaotic properties probably including expo- 
nential decay of multiple correlations [26,27]. In the case 
R — 1/2 the circles touch tangentially leading to the possi- 
bility of long sequences of collisions near the corners (hence 
weaker chaotic properties) , and the limit R — ► 1/ \[2 is the 
intcgrable (regular) square. Here we consider R = l/y/3 
which leads to simple exact formulas for the perimeter 

|d£>| = 2tt/3\/3, area \D\ = 1 - ^4/^3-1 - tt/9 and 
mean free path (T) = Tgjjl- The latter formula is general 
for two dimensional billiards [28] . The hole position coor- 
dinate I* is defined so that a curved side corresponds to a 
unit interval. 

For the numerical simulation, a trajectory of 10 8 colli- 
sions is simulated and stored. The position on the bound- 
ary is binned according to a partition of the arcs into pieces 
of fixed length, giving a single sequence of integers, from 
which escape rates of open systems with holes given by 
any desired combinations of the pieces can be calculated 
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Fig. 2: (Color online; monochrome versions have from light 
to dark: yellow, green, red, blue) The escape rate 7 (red), as 
a function of the position of two holes of size R/20, together 
with the first (blue) and second (green) order approximations 
in Eq. 19 and using one-hole escape rates (yellow) as in Eq. 21. 
Note that the 3D graphics format does not conceal large fluc- 
tuations, for example the RMS deviations from the directly 
calculated escape rate are 0.00848 (first order), 0.00193 (sec- 
ond order) and 0.00190 (second order using one-hole escape 
rates) . 
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Fig. 3: The escape rate 7 (solid) as a function of the position 
of one hole of size R/50 at position I*, together with the first 
(dotted) and second (dashed) order approximations in Eq. 19. 

in regular systems [19]. 

L.B. was partially supported by NSF grant DMS- 
0140165. 



simultaneously. The correlation functions are also calcu- 
lated from the time series; for the infinite sums only a 
few (roughly ten) terms need to be retained, due to the 
exponential decay of correlations for this system. 

The choice of the size of the hole affects the numerical 
tests in that for large holes higher order (h 3 ) terms be- 
come significant, while for small holes the h 2 terms may 
be smaller than the errors in the correlation statistics. 

The results of the numerical simulations for some hole 
sizes and positions are shown in Figs 2-4. Clearly the sec- 
ond order approximations are better than the first order 
approximation, especially when both holes are near the 
same corner, which leads to a lower escape rate. Near 
a corner, we expect that some orbits have several colli- 
sions in the holes, and hence the higher order correlation 
functions are more important. In Fig. 3 we note that 
fluctuations in the escape rate as a function of hole 
position are visible; compare with Ref. [15]. Fig. 4 
gives useful information about the validity of the 
expansion truncated at the first or second term 
even when the hole is no longer small. 

The numerical simulation thus confirm the formulas to 
second order, which is already sufficient to determine de- 
tailed information about the internal dynamics from ob- 
serving escape rates as a function of hole position. Com- 
putation of the necessary higher cumulants where holes 
are on short periodic orbits will require a deeper analysis 
of the role of the latter. Another important but only par- 
tially understood problem is that of general systems with 
power law escape, such as stadium-type chaotic billiards, 
and systems with mixed dynamics, where chaotic and reg- 
ular behavior coexist. We are confident that our approach 
(suitably modified) will work and lead to general and use- 
ful formulas as well, as suggested by similarities between 
results in the strongly chaotic systems considered here and 
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Fig. 4: As for Fig. 3, except that the hole is of variable size, 
extending from zero to I* . 
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